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Abstract 



We present spherically symmetric black hole solutions for Einstein gravity 
coupled to anisotropic matter. We show that these black holes have arbitrar- 
ily short hair, and argue for stability by showing that they can arise from 
dynamical collapse. We also show that a recent 'no short hair' theorem does 
not apply to these solutions. 
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I. INTRODUCTION 



It is widely believed that black holes retain only very limited information about the mat- 
ter that collapsed to form them. This information is reflected in the number of parameters 
characterizing the black hole solution. For stationary axis-symmetric black holes in general 
relativity with no matter coupling, the only such parameters are the mass M and angu- 
lar momentum J. With coupling to electromagnetism, the parameters include the electric 
charge Q. These results are known collectively as the 'no hair' theorems JO-PI. 

The three parameters M, J, and Q are special in that they can be detected by performing 
experiments in the asymptotically flat region of the black hole spacetime. This is because 
they are captured by conserved surface integrals in the asymptotically flat region. They are 
not what is referred to as 'hair'. Any parameters in a black hole spacetime that are not 
captured by asymptotic surface integrals may be called 'hair'. Thus, no hair theorems imply 
that all the parameters characterizing black hole spacetimes are those that are 'visible' at 
spatial infinity. 

For Einstein gravity coupled to matter fields other than electromagetism, black holes 
might be parametrized by quantities other than M, J, and Q. In these cases, the parameters 
might all be captured by surface integrals at infinity, leading to an associated 'no hair 
theorem'. Alternatively, such theories might yield black hole solutions with hair. 

In this paper we describe spherically symmetric and static black hole solutions that arise 
from coupling Einstein gravity to an 'anisotropic fluid', which is defined below. For certain 
special cases this matter gives the de-Sitter, and Reissner-Nordstrom black holes. For other 
cases, the resulting black holes have hair in the sense described above. In the next section 
we describe the anisotropic medium and construct a class of static, spherically summetric 
solutions to the coupled Einstein-matter equations of motion. In section III we specialize to 
a class of black hole solutions and show that these black holes have short hair. We also give 
arguments indicating that these black holes are stable. Section IV discusses the relationship 
of our black hole solutions with the so called 'no short hair' conjecture of Ref. 0. We 
conclude in section V with a short discussion. 



II. SOLUTIONS WITH ANISOTROPIC MATTER 

Consider a static, spherically symmetric metric of the form 

ds 2 = -f(r)dt 2 + -^--dr 2 + r 2 d6 2 + r 2 sin 2 6d(j) 2 (2.1) 
f( r ) 

The nonvanishing components of the Einstein tensor for this metric are 

rf'-l + f 



r 



Lr t - U r - — 2 

G e e = G+ = , (2.2) 

where the prime denotes a derivative with respect to r. It follows that any matter source 
that gives rise to a metric of the form Eq. (2.1) must satisfy T\ = T£ and T$ = Tf. 
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In particular, the "radial pressure" T r r need not equal the "angular pressure" Tff = 
Therefore, generically, such a matter source is not isotropic in the static time slices. What 
we seek, then, is a medium that is capable of supporting anisotropic stresses. 

A material with the desired properties is the elastic medium discussed by DeWitt [|J. 
We begin with DeWitt's action, as written in Ref. |§, 

S[g ab , Z l \ = - f d A x^g~ p(Z\ h jk ) . (2.3) 

This action is a functional of the spacetime metric g ab and the Lagrangian coordinates Z l (x), 
i = 1,2, 3. The Lagrangian coordinates are a set of labels C = Z l {x) that tell which particle 
passes through the spacetime point x. Thus, Z % (x) are three scalar fields whose gradients 
are orthogonal to the matter world lines. The action S[g a b, Z % \ is the proper volume integral 
of the Lagrangian — p, where p is the proper energy density in the rest frame of the matter. 
The energy density p depends on Z l explicitly, and also implicitly through the matter space 
metric h^. The matter space metric is defined by 

h* = {d a V)g a \d h Zi) , (2.4) 

and is interpreted as the metric in the rest frame of the matter. That is, hijdQdC? is the 
square of the proper orthogonal distance separating particle world lines with Lagrangian 
coordinates Q and Q + dQ . Using hij, the spacetime metric g a b may be written in the form 

g ab = -u a u b + h t] {d a z l ){d b z3) , (2.5) 

where U a is the fluid velocity, defined as the future pointing unit vector orthogonal to the 
gradients d a Z\ 

The stress-energy-momentum tensor for the elastic medium is 

r "> s = ~ pg - + 2 ^ ld " z ' mz,) 

= pUJJ b + t, J (daZ i )(d b Z i ) , (2.6) 
where we have used (12.51) and defined the matter stress tensor 



dp \ 2 d(Vhp ) 



= 2 t£; ~ A, ] = — • (2-7) 



The matter equations of motion, obtained by varying the action with respect to Z l , are 
equivalent to the conservation equations V a T a6 = 0. When the Einstein equations hold, the 
contracted Bianchi identity implies that the matter equations are also satisfied. 

An isotropic elastic medium is one whose principal pressures are all equal. Its stress 
tensor is that of a perfect fluid, and is obtained by letting p depend only on the proper 
volume of the Lagrangian coordinate cell d 3 ( and the number of particles in that cell. This 
is realized by setting 

p = p(N) , (2.8) 
with the scalar iV defined by the ratio of densities 
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N = N{Z l )/Vh , (2.9) 

where N_(C) is the number of particles per unit coordinate cell d 3 (. Thus, N is the proper 
number density, i.e., the number of particles per unit proper volume in the material rest 
frame. With these definitions, the stress-energy-momentum tensor (|2.6| ) assumes the famil- 
iar form: Using the identity 

dfri 2 ON 11 ' 1 UJ 

and (|2.5|), the tensor ( |2.6| ) becomes 



T ab = P U a U b + - p^j (g ab + U a U b ) . (2.11) 

Thus, the matter action ( |2.3|) with p = p(N) is the action of a perfect fluid with pressure 
P = N(dp/dN) — p. This definition of pressure can be written as d(p/N) = —Pd(l/N), 
which is the first law of thermodynamics applied to a cell with a fixed number of particles. 

An anisotropic elastic medium is one whose principal pressures are not all equal. If two 
of the three principal pressures are equal, the fluid has planar symmetry. The basic idea 
for introducing anisotropy is to perform a '2 + 1' decomposition of the matter space metric 
hij, and treat the plane differently from the line. This may be done as follows. Let us view 
the Lagrangian coordinates ( l as points in a manifold S, which is the space of particle flow 
lines, and introduce a foliation of S by a function <^(0) with d{ip ^ 0. Now set p = p(n, m), 
where n is the particle number per unit proper area in the ip = const slices and m is the 
particle number per unit proper length in the direction orthogonal to these slices. This 
introduces the anisotropy. For definiteness and simplicity, let us choose <p(C) — C 1 - Then 
n = n/ \fa where n(0) is the number of particles per unit coordinate area d 2 ( = dQ 2 A d( 3 , 
and a = ^22^33 — (^23) 2 is the determinant of the metric in the C 1 = const slices. Similarly, 
m = m/a where m is the number of particles per unit coordinate length d^ 1 and a = Jh/cr 
is the 'lapse function' that measures proper length in the direction orthogonal to Q 1 = const. 

We thus have the following '2 + 1' form for h^: 

where (3^ is the 'shift vector' and the indices /1, v take the values 1, 2. Equivalently, we have 
= UiUj + o~ij, where U{ = ab~} is the unit vector orthogonal to the ( l = const, surfaces. 
The corresponding spacetime tensors are 

U a = U^aZ 1 = adaZ 1 , CJ ab = (hij - u i u j )(d a Z i )(d b Z') , (2.13) 

which represent, respectively, the unit normal and metric of the isotropic 2-surfaces in the 
rest frame of the matter. 

The stress tensor for the anisotropic medium, with p(n,m), is obtained by computing 
(dp/dh l i) using the above definitions, and substituting into ( |2.7| ). We find 

dp(n,m) ndp m dp . _ , s 

-M = 2 ^" + (2 - M) 
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which gives 



Finally, the full stress-energy-momentum tensor for the anisotropic matter is 

Tab = pU a U b + ~ P^j °ab + ~ U a U b • ( 2 -16) 

We note, as required, that the pressure in the 2-surfaces is n(dp/dn) — p, which differs from 
the pressure m(dp/dm) — p in the direction orthogonal to 2-surfaces. 

Our goal now is to find an equation of state, p = p(n,m), that will yield the metric 
(27T) as a solution of the Einstein and matter equations of motion, for arbitrary f(r). The 
Einstein tensor ( |2.2|) shows that the spheres t = const, r = const are isotropic. Thus, we 
seek a solution where u a is radial. In particular, let Z x {x) = r, Z 2 (x) = 8, and Z 3 (x) = <fi 
so that u a = a5 r a . From ( |2.16|) the condition T\ = T r r now shows that the equation of state 
must be independent of m, so p = p{n). In this case we have 

Ti = Tl = -p (2.17) 

T °= T t = n ^-P> ( 2 - 18 ) 
and the Einstein equations yield (with G = 1) 

rf'-l + f 



r 



2 



up (2.19) 



rf" + 2f (dp \ 

-^— = M \ n d^- p ) ■ (2 ' 20) 

Using the relationship {rf — 1 + /)' = rf" + 2/' between the numerators of the left-hand 
sides of these equations, we find 

(—r 2 p)' dp 

n^--p. (2.21) 



2r dn 
This reduces to n'/n = — 2/r, and implies 

n = £ , (2.22) 

where K is a positive dimensionless constant. Now define the function F(n) by 

F(n) = F(K/r 2 ) = f(r) . (2.23) 
Then the Einstein equation ( 2.19 ) gives 

n ( dF\ 
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This determines the equation of state p(n) in terms of the arbitrary function f(r) that 



appears in the metric ( 2.1 



Note that the correspondence between f(r) and p(n) is not one-to-one. One can see 
this by solving Eq. Q2.24Q for F(n): 

F{n) = 1 + ^Kyfn J dnn- 5/2 p(n) . (2.25) 

The integral over n yields a term AirKK' y/n which is proportional to the integration constant 
K' . This term corresponds to — 2M/r in the metric function /(r), where M = —2tiK 3 / 2 K > . 



Let us reexpress and summarize these results. Consider an elastic medium fl2.3H2.4|) with 
equation of state p = p(n), where n = n/^/a, n = K sin Z 2 , and a = ^22^33 ~ (^23) 2 - A 
family of solutions of the Einstein and matter equations of motion is Z 1 (x) = r, Z 2 (x) = 6, 
Z 3 (x) = 4>, and g a b(x), where g a b(x) is given by (|2.1| ). The function /(r) that appears in 
the metric is defined by f(r) = F(K/r 2 ), where F{n) is determined by the equation of state 



through Eq. ( |2.25|) . This family of solutions is parametrized by an arbitrary constant M 



that appears in a term — 2M/r in the metric function f(r). Equation ( |2.24| ) shows that, 
for a suitable equation of state p = p(n), any metric of the form (|2.1| ) can be obtained as a 
solution of the Einstein and matter equations of motion. 

The anisotropic elastic medium described here satisfies the (weak, dominant, or strong) 
energy conditions only if restrictions are placed on the equation of state. For example, 
for the weak energy condition the proper energy density must be nonnegative, p > 0, and 
the principle stresses must be greater than or equal to —p. This later requirement implies 
dp/ On > 0. The inequalities p > and dp/ On > can be translated into restrictions on the 
function F(n), which in turn restrict the form of the metric function f(r). 



III. BLACK HOLES WITH SHORT HAIR 

In the previous section we showed that the general spherically symmetric metric ( [2.1| ) 
arises from an anisotropic matter source. However all such metrics are not physical because, 
as mentioned above, the matter does not in general satisfy physical energy conditions. In 
this section we consider a class of black hole solutions of the Einstein-anisotropic fluid theory 
that are physical. We show that these black holes have arbitrarily short hair, and argue that 
they are stable under scalar perturbations. The stability argument relies on the observation 
that these black holes can arise from spherically symmetric gravitational collapse. 

Consider the equation of state 

p(n) = Cn k+1 (3.1) 



for the anisotropic fluid, where C and k are constants. From Eq. (|2.16 ), the stress-energy- 
momentum tensor is 

Tab = pUJJb - pu a u b + Pa ab , (3.2) 

where 

P = n^-p = kp (3.3) 
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defines the principal pressures in the 6 — <p surfaces. 

T ab is a type I stress tensor, so the weak energy condition requires p > and P > —p. 
This clearly holds for all C > 0, k > — 1. The strong energy condition requires, in addition to 
the weak energy condition, that the sum of the energy density p and the principal pressures 
should be non-negative. In our case this implies P > 0, so the strong energy condition 
holds for all C > 0, k > 0. The dominant energy condition requires, in addition to the 
weak energy condition, that energy fluxes should never be spacelike. In our case this leads 
to p > P. Therefore the dominant energy condition holds for C > 0, — 1 < k < 1. 

From Eq. ( |2.25| ) we see that the metric resulting from the equation of state p(n) = Cn k+l 
is O) with 

/ 2M 2k \ 

ft) ={ 1 - — + 7^) ' ^ ^ ■ ^ 

Here, we have chosen C > and defined Q 2k = 8irCK k+1 /(2k — 1). The energy density and 
pressures are 

p _ Q 2k (2k-l) p _ h Q 2k (2k-l) 

Notice that for k — 1 this solution is identical to the Reissner-Nordstrom metric. For our 
purposes, the more interesting situation occurs for k > 1. In that case the parameter Q 
is not captured by a surface integral at spatial infinity due to the rapid fall off r~ 2k of 
the corresponding term in f(r). Therefore, when k > 1, Q qualifies as hair. Since we 
can take k to be arbitrarily large, the hair Q can be arbitrarily short. This means that 
the dependence of the metric on Q may be confined to an arbitrarily small neighbourhood 
of the curvature singularity at r = 0. As an example, consider the limit Q <C M. In 
this case the outer event horizon is at r + pa 2M and the inner (Cauchy) horizon is at 
r- ~ <3(<3/2M) 1/(2fc-1) + 2M(Q/2M) 4k ^ 2k -^/(2k - 1). The ratio of the energy densities at 
the event and inner horizons is p(r + )/ p(r_) = (r_/r + ) 2fc+2 w (Q/2M) 2fc ( 2fc+2 )/( 2fc_1 \ This 
ratio can be made arbitrarily small by decreasing Q/2M or by increasing k. Thus the hair 
sprouting from the event horizon can be arbitrarily short. We emphasize that the anisotropic 
fluid that gives rise to these short-hair solutions satisfies both the weak and strong energy 
conditions. 

We now consider the issue of stability for the black holes with short hair. In particular, 
we present evidence that the black holes ( |3.4| ) are stable under the spherically symmetric 
perturbations 

/(r)-/(r) + <f/(r,t), p(r) - p(r) + 5p(r,t) , P(r) -> P(r) + 6P(r,t) (3.6) 

about the static solutions. As discussed below, these are not the most general possible 
spherically symmetric perturbations. 

Stability of a static or stationary black hole is studied by performing a standard but te- 
dious analysis to obtain the equations governing metric and fluid perturbations to linearized 
order. The black hole is stable if initial metric perturbations do not grow without bound. 
Alternatively, one can attempt to show that the black hole arises from the long time limit 
of gravitational collapse. Indeed, if a static or stationary black hole is known to arise as the 
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endpoint of a dynamical evolution, its stability is not seriously in question. We will consider 
an argument for stability based on the second approach; however, it should be emphasized 
that such an approach is physical, and is not a substitute for a complete mathematical 
analysis of the first type. 

To begin, let us observe that the stress-energy-momentum tensor ( |3. 2| ) can be rewritten 

as 



Tab = p(r)(v a w b + v b w a ) + P(r)(g ab + v a w b + v b w a ) , (3.7) 
where the null vectors v a and w a are defined by 

U a = ~^(v a + w a ) , U a = —f=(v a - W a ) . (3.8) 

For our solutions U a = (^,—1/^,0,0) and u a = (0, 1/a/J, 0, 0) in the coordinates 
(v, r, 9, 0), where v is an advanced time coordinate defined by dv = dt + (1/ f)dr. 

In [TO] a stress tensor similar to ( |3.7| ) was used to obtain dynamical spherically symmetric 
collapse solutions. That tensor is 

Tab = , \ v a v b + p(r, v) (v a w b + v b w a ) + P(r, v) (g ab + v a w b + v b w a ) . (3.9) 

The general solution of the Einstein equations for this source, with the equation of state 
P = kp, is 

dS* = -(l- + 2 9(V) \ ^2 + 2dvdr + (310) 



where 



r (2k — l)r 



m(r,v) = f(v) - {2k ^ 2k _ t , (Ml/2), (3.11) 



nr,v) = k£^- 2 = kp(r,v), (3.12) 

and f(v) and g(v) are arbitrary functions. Now we can choose the functions / and g such 
that lim^oo f(v) = M and lim^ 00 ^(w) = Q 2k (2k - l)/2. Then the metric ( pTTCD at large 
values of advanced time coincides with the short hair black hole Q2.1] ), ( |3.4|) . In this way we 
see that the static short hair solutions result as the end point of gravitational collapse. The 



metric ( |3.10|) is to the short hair black hole what the Vaidya metric is to the Schwarzschild 
black hole. 

We now elaborate on the connection between the stress-energy-momentum tensors ( |3.7|) 
and ( |3.9| ) induced by the perturbations (|3.6|). Notice first that if the functions f,p and P 
are all taken to have general (r, v) dependence, the most general stress-energy-momentum 
tensor must be of the form ( |3.9| ) JTD|. Thus the perturbations (p.6[), when substituted into 
the tensor (|3.7|), must lead to a tensor of the general form (|3.9|); the extra v a v b term is 



effectively induced by the time dependence of /. Now since we have the general solution 
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( ft.lOp with the source ( |3.9| ), we have implicitly the solution of the equations governing the 
perturbations ( |3.6|) , to all orders in the perturbations. This is because the general solution 
( p.lO| ) can always be separated into a static part plus a time dependent 'perturbation'. Thus, 
this argument shows stability under perturbations of the type (|3.6|) . 

The most general spherical scalar perturbations are however not given by ( |3.6| ). To 
obtain the most general perturbations of this type we would need to add to the stress- 
energy-momentum tensor ( |3.9| ) an 'outgoing' (r, v) dependent term proportional to w a Wb. 
Such a term comes effectively from including a metric perturbation 8h(r, v) dvdr. Although 
an exact collapse solution with some matter outgoing to infinity is not known, one can 
argue, by analogy with the numerical spherically symmetric scalar field collapse fllTf , that 
an outgoing component with positive energy density will not affect stability when black 
hole formation occurs. In that system it is observed, in the so called 'supercritical' region 
of parameter space (where a black hole forms), that an initially ingoing pulse oscillates 
radially as it collapses, with some scalar waves escaping to infinity with each oscillation as 
the collapse to a black hole proceeds. Similar behaviour has been observed numerically in the 
collapse of perfect fluids |r2| , and may occur for the type of anisotropic fluid discussed here. 
While this analogy is suggestive, we only have an argument for stability for the spherical 
'ingoing' perturbations 



IV. COMMENT ON A 'NO SHORT HAIR' THEOREM 

As we have seen, some of the static spherically symmetric black hole solutions (those 
for k > 1) arising from gravitational coupling to anisotropic media have hair — that is, 
parameters in the metric which are not captured by surface integrals at spatial infinity. The 
energy density falls to zero with radius r faster than 1/r 4 . The physical picture is that some 
matter protrudes from the black hole event horizon without being pulled in, but does not 
'extend to infinity' as in the electric field case. Matter can hover in a strong gravitational 
field without collapsing completely only if its internal pressures are sufficiently large. This 
appears to be the case for the anisotropic fluid black holes with k > 1. 

More generally one can ask: How close to a black hole event horizon can matter hover 
if the stress-energy-momentum tensor satisfies some reasonable physical conditions? Or in 
other words, how 'short' can hair be? This question has been posed, and there is a theorem 
based on certain assumptions |7| . Specifically, the theorem states that if matter is such that 

(i) the weak energy condition holds, 

(ii) the energy density p falls to zero faster than r~ 4 , 

(iii) T := T ab g ab < 0, 

then it follows that Pr A (P = radial pressure) is negative and decreasing until at least 3/2 
the radius of the event horizon, beyond which this quantity increases to zero. This behaviour 
is interpreted as 'long hair', and it is suggested ]7| that it will occur for any non-linear 
theory coupled to gravity. The three conditions above effectively imply that radial pressure 
is negative, at least near the horizon. 

We have presented black hole solutions that have arbitrarily short hair for the parameter 
k > 1. For our 'short hair' black holes, conditions (i) and (ii) hold, but condition (iii) does 
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not: T = 2p[k — 1), which is positive for k > 1. Thus, our solutions do not contradict the 
no short hair theorem because the condition T < for it to hold is not satisfied. Thus, 
while this theorem does apply to several discussed in 0, it does not apply to the 

anisotropic fluid. 



V. DISCUSSION 

We have presented a class of black hole solutions arising from Einstein gravity coupled 
to a type of anisotropic matter, and described some of their properties. Perhaps their most 
unusual feature is the possibility of arbitrarily short hair. 

It is interesting to observe that these black holes can also be seen as arising from a 
generalized electromagnetism. Consider the stress-energy-momentum tensor 

T ab = i (F ac F b c - | g ab F cd F cd ) , (5.1) 

where F ab is the field strength and a is a constant The field equation for F ab is that 
obtained from V a T afe = 0. A solution of the Einstein equations coupled to this matter is 
the metric Q2.1|) and (|3.4[) , with 



ab 



' (a-iXE + i) ^ {dtAdrUt (5 . 2) 



and a = 2k/ (k + 1). We note the black hole solutions viewed as arising from this matter 
source can be shown to be stable via an analysis very similar to that given in |H| for 
the Reissner-Nordstrom solution. Indeed, the only difference in the details are appropriate 
replacements of Q 2 /r 2 by Q 2k /r 2k . 

It would be of interest to study the collapse of anisotropic fluids numerically, especially 
in light of the black hole solutions we have found. In particular, as has been done for the 



perfect fluid ||12|| , one could seek a 'critical' solution with a self-similarity ansatz. If found, 
such a solution could be used as the starting point of perturbation theory pl2] , |T4|| for finding 
critical exponents associated with the parameters M and Q. 

Note added: After completing this work we became aware of related results by Magli and 



Kijowski [15|. These authors consider equilibriuim, spherically symmetric stars composed 



of anisotropic elastic matter. We thank G. Magli for bringing this to our attention. 
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We do not know an action which leads to the stress-energy-momentum tensor fl5.1|) . 
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